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ABSTRACT 

Aims. We present a survey of the ortho-H2D^(li o-li,i) line toward a sample of 10 starless cores 
and 6 protostellar cores, carried out at the Caltech Submillimeter Observatory. The high diag- 
nostic power of this line is revealed for the study of the chemistry, and the evolutionary and 
dynamical status of low-mass dense cores. 

Methods. The derived ortho-H^D^ column densities (N{ortho - HaD^)) are compared with pre- 
dictions from simple chemical models of centrally concentrated cloud cores. 
Results. The line is detected in 7 starless cores and in 4 protostellar cores. N{ortho - H2D*) 
ranges between 2 and 40x10'^ cm"^ in starless cores and between 2 and 9x10'^ cm"^ in proto- 
stellar cores. The brightest lines are detected toward the densest and most centrally concentrated 
starless cores, where the CO depletion factor and the deuterium fractionation are also largest. The 
large scatter observed in plots of N{ortho - HiD^) vs. the observed deuterium fractionation and 
vs. the CO depletion factor is likely to be due to variations in the ortho-to-para (o/p) ratio of 
H2D+ from > 0.5 for < 10 K gas in pre-stellar cores to - 0.03 (consistent with T\^\^ ^ 15 K 
for protostellar cores). The two Ophiuchus cores in our sample also require a relatively low o/p 
ratio 0.3). Other parameters, including the cosmic-ray ionization rate, the CO depletion factor 
(or, more in general, the depletion factor of neutral species), the volume density, the fraction of 
dust grains and PAHs also largely affect the ortho-H2D^ abundance. In particular, gas temper- 
atures above 15 K, low CO depletion factors and large abundance of negatively charged small 
dust grains or PAHs drastically reduce the deuterium fractionations to values inconsistent with 
those observed toward pre-stellar and protostellar cores. The most deuterated and H2D*-rich 
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objects (L 429, L 1544, L 694-2 and L 183) are reproduced by chemical models of centrally con- 
centrated (central densties ^10* cm"^^) cores with chemical ages between 10* and 10* yr. Upper 
limits of the para-H3 0^(l7-2|) and para~D2H^(l i.o-lo,i) lines are also given. The upper limit to 
the para-HjO^ fractional abundance is ^ 10 ** and we find an upper limit to the para-D2H+/ortho- 
H2D^ column density ratio equal to 1, consistent with chemical model predictions of high density 
(2x10* cm"^) and low temperature (Tii„ < 10 K) clouds. 

Conclusions. Our results point out the need for better determinations of temperature and density 
profiles in dense cores as well as for observations of para-H2D^. 

Key words, astrochemistry — stars: formation — ISM: clouds, molecules — radio lines: ISM 
— submillimeter 



1. Introduction 



In the past decade, astrochemistry has become more and more crucial in understanding the structure 
and evolution of star forming regions. There are no doubts that stars like our Sun form in gas 
and dust condensations within molecular clouds and that the process of star formation can only 
be understood by means of detailed observations of the dust (good probe of the most abundant 
and elusive molecule, H2) and molecular lines (unique tools to study kinematics and the chemical 
composition). 

Millimeter and subm illimeter continuum dust emission observations (see 
Ward-Thompson et al. 1 12007[ for a detailed review of this topic) have provided a very good 



probe of the density structure of dense cores, although uncertainties are still present regarding 
the dust opacity and temperature, both likely to change within centrally conce ntrated objects 
(but such variations have so fa r been hard to quantify observationally, see e.g. 



Bianchi et al. 



2003 



Pagani et al. 



2003 



20041) . Stellar counts in the near - infrared provide an alternative way of 



measuring the dust (and H2) column and cloud structure (ILada et al. Ill994l) . independent of any 
variat ion in dust proper ties, but they cannot probe regions with extinctions above about 40-50 
mag (lAlves et al. Ill998), i.e. the central zone of ver y dense cores, such as L 1544, where Ay - 
100 mag within 11 jWard-Thompson et al. Ill999h. It appears that many starless cores can be 



1955 



Bonnor 



approximated as Bonnor-Ebert spheres dEbert 
densities ranging from about 10^ cm~^ (as in the case of B 68; 



e.g., L 1544, L 183, and L 694-2: IWard-Thompson et al. 



1956h. with valu es of the central 
Alves etal. 11200 lb to lO '' cm'^ (for 



1999 



Pagani et al. 



2003 



Harvev et al. 



2003b . At the lower end of t he central density range, dense cores ap pear to be isothermal, with gas 



temperatures close to 10 K (IGalli et al. 



2002 



Tafalla et al. 



20041) . whereas higher density cores 



have clear evidences of temperature drops in the central few thousand AU, with dust temperatures 



approaching about 7 K ( 


Evans et al. 


2001 


Pasani et al. 


2003 


2004 


Schnee & Goodman 


2005 


Pagani et al. 


2007 


Crapsi et al. 


200' 


, see also Bergin & Tafalla 2007 for a comprehensive review 



on starless cores) 



Keto & Field 



(120051) have proposed that the "shallower" cores (such as B 68) are in approxi- 
mate equilibrium and will not evolve to form protostars, whereas the centrally concentrated ones 
(such as L 1544) are unstable cores that are proceeding toward gravitati onal collapse and th e for- 
mation of protostars. Indeed, this is in agreement with the findings of 
claim that B 68 is oscillating around an equilibrium state, and those of 



Lada et al. 



Caselli et al. 



(12002). who 



(l2002al) and 
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van der Tak et al. 



(I2005h . who studied the kinematic structure of L 1544 and found that it is con- 
sistent with contraction in the core nucleus (or central contraction). 

To trace the gas pro perties, NH3 and N2H^ have been extensively used for several years 



dBenson & Mvers 



19891) an d they seem to trace quite similar conditions, having comparable mor 



phologies and line widths (IBenson et al. 



1998 



Caselli et al. 



2002c : 



Tafalla et al. 



2002, 



2004 . 



despite of the (two orders of magnitude) difference i n the critical densiti es of the most frequently 
observed transitions (NH3(1,1) and N2H+(l-0), see IPagani et al. 1 120071) . This is quite lucky for 



astronomers considering that only relatively recently it has been realized that these two species 
are among the few that are left in the gas phase at volume densities above ~10^ cm"^. In fact, 
CO, C S and, in general, all the carbon bearing species so far observed (with the exception 



of CN; 



Hilv-Blant et al 



starless cores 



Kuiper et al. 



20081) are heavily affected by freeze-out in the central parts of dense 



1996; 



WiUacv et al. 



1998 



CaseUi et al. 



1999; 



Herein et al. 



2001 



Bacmann et al. 12002 ; 



Herein et al. 12002: 



Caselh et al. l2002b 



Tafalla etal. 12002: 



Hacmann et al. 



2003 



2006; 



1984 



Crapsi et al. 



2004 : lTafalla et al. ll2004lPagani et al. IbOOSHCrapsi et al. ll2005l:lTafalla et al 



Pagani et al. 



Roberts & Millar 



2007 ). The f reeze-out of neutral species (in particular of CO. 



200o3bl) . boosts the deuterium fractionation in species such as N2H^ , 



Dalgarno & Lepp 



NH3, H7CO , and HCO^(see e.g 



Hacmann et al. 



2003 



Crapsi et al. 



Hutner et al. 



2005 



1995 



Lis et al. 



Tine et al. 



2006 



200c ; 



Gerin et al. 



Caselli et al. 



2002b 



2006). Also in star form- 



ing cores, in part icular in the direction of Class sources, the first protostellar stage (e.g. 



Andre et al. 


2000 


, right after the pre-stellar phase, the deuterium fractionation is found to be 


verv large ( Ceccarelh et al. Il998 


Loinard et al. I2OO2 


:Lis etal.l2002jb; 


van der Tak et al. II2OO2; 


Parise et al. 


2002 


Vastel et al. 


2003: I 


'arise et al. 


2004: Crapsi etal. 


2004; Marcelino et al. 


2005: Parise etal. 


2006: Ceccarelli et al 


2007). This is thought to be the signature of a (recent) 



past in which the star forming cloud core experienced the low temperature and high density con- 
ditions typical of the most centrally concentrated starless cores, where some deuterated molecules 
(e.g. deuterated ammonia) are formed in the gas phase (and stored on dust surfaces) whereas oth- 
ers (such as deuterated methanol and formaldehyde) are likely formed onto dus t surfaces and then 



released to the gas phase partially upon formation dOarrod et al. 



2006 . 



2007h as well as via the 



Turner 



199C; 



Cazaux et al. 112003 



interaction with the newly born proto star, which can (i) heat dust grains, leading to mantle eva .p- 
oration (as in Hot Cores and Corinos; 



Hottinelh et al. Il2004allb 



I 



20071), and/or (ii) " erode" dust man tles via sputtering in shocks produced by the associated ener- 



getic outflows (e.g. 



. Lis et al. Il2002al) . 

Questions that are still open are: (1) are N2H^, NH3 and their deuterated forms really tracing 
the inner portions of centrally concentrated cores on the verge of star formation? Although Hergin 
et al. (2002), Pagani et al. (2005, 2007) found evidences of depletion in the center of H 68 and 
L 183, respectively, there are no signs of freeze-out for NH3 in L 1544 (Crapsi et al. 2007) and 
for both nitrogen bearing species in L 1517B and L 1498 (Tafalla et al. 2004). At densities above 
~10^ cm"^, the freeze-out time scale is quite short (~1,000 yr) and all heavy species are expected 
to condense onto grain mantles. Moreover, recent laboratory measurements clearly show that N2, 
the parent species of both NH3 and N2H^, s hould freeze-out at the same rate as CO (h aving similar 



binding energies and sticking probabilities: lOberg et al. 



2005 



Hisschop et al. 1120061) . (2) For how 



long the high degree of deuterium fractionation observed in starless cores is maintained after the 
formation of a protostellar object? 
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The detection of strong ortho-H2D^(l 1,0 - 1 1,1) emission in the direction of L 1544 



dCaselU et al. 



20031) . and the conclusion that Ht, with its deuterated counterparts, is one of the 



most abun dant molecular ions in co r e centers, have opened a new way to study the chemic al 



evolution ( 


Roberts et al. 


2003 


2004; 


Waknslev et al. 


2004 


; Flower et al. 


2004, 


2005, 


2006*: 


Aikawa et al. 


2005 


) and the kinematics ( 


van der Tak et al. ^ 


1005 


) of the central few thousand AU 



of starless cores. Thus, H2D^ is an important tool to understand the chemical and physical proper- 
ties of the material out of which protoplanetary disks and ultimately planetary systems form. 

In the gas phase at temperatures below ~20 K, the deuterium fractionation is mostly regulated 
by the proton-deuteron exchange reaction: 



H+ + HD ^ HzD^ + H2 + AE, 



(1) 



where AE 230 K, 



Millar et al. 



19891) prevents th e reverse reaction to be fast in cold regions 



unless a sign i ficant fr action of H2 is in ortho form dGerlich et al. 



2002; 



Walmslev et al. 



Flower et al. 



2004) 



(l2006al) showed that values of the ortho-to-para (o/p) H2 ratio much higher than 
0.03 are inconsistent with the observed high levels of deuteration of the gas (see their Fig. 6). The 
above reaction, together with the fr eeze-out of neutral species (which boosts up the production 
rate of H2D^ compared to Ht; e.g. 



Aikawa et al. 



200 Ih . allows the H2D^/Ho ratio to overcome 



the cosmic D/H ratio by several orders of magnitude. In fact, not only ortho-H2D^(li_o - li,i) has 
been found to be ^ I K strong in L 154 4, but also D2H^ ha s been detected toward another starless 
condensation (16293E, in Ophiuchus; Vastel et al. 2004 1. The ortho-H2D^(l 1 - li,i) line has 



also recently been mapped in L 1544 (iVastel et al. Il2006al) . finding that the H2D^ emitting region 
has a rad ius of about 5000 AU , comparable to the size of the N2D^(2-1) map made in the same 
region by 



Caselli et al. 



(I2002ah . 

In this paper we present new ortho-H2D^(l 1 - li,i) observations, carried out at the Caltech 
Submillimeter Observatory (CSO) antenna, in the direction of 10 starless cores and 6 cores asso- 
ciated with very young protostellar objects. As it will be shown, the line has been detected in 7 of 
the 10 starless cores and in 4 out of 6 star forming cores. In Sect. |2]the observational details are 
given. Results and ortho-H2D^ spectra are shown in Sect. [3] together with a brief discussion on the 
upper limits of the para-D2H^(l i o-lo.i) lines (para-H30^(l j^-2|) upper limits can be found in the 
on-line appendix). H2D^ column densities are derived in Sect.|4] A chemical discussion, aimed at 
interpreting the observations, is described in Sect.|5]and conclusions are in Sect.|6] 



2. Observations 

2.1. Technical details 

Obser vations of the ortho-H2D+(li,o - li,i) line (vq = 372.421385(10) GHz; 



Amano & Hirao 



20051) were carried out at the Caltech Submillimeter Observatory (CSO) on Mauna Kea (Hawaii), 
between October 2002 and April 2005. The spectra were taken in wobbler switching mode, with a 
chop throw of 300". The backend used was an acousto-optical spectrometer (AOS) with 50 MHz 
bandwidth. The velocity resolution, as measured from a frequency comb scan, is 0.1 kms The 
beam efficiency (rib) at v = 372 GHz was measured on Saturn, Mars and Jupiter and is listed 
in Tab. |2] Measurements for extended sources were made for only few sources (L 1544, L 183, 
NGC 1333-DCO+, B 1, NGC 2264G-VLA2) and found to be ~ 70%, compared to 60% on planets 
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measur ements. HaD^is likely to be extended (e.g. L 1544: 



Vastel et al. 


2006a 


;L 183: 


Vastel et al. 



2006bl) . Consequently, we used the extended source beam efficiency whenever available. However 



note that the difference cannot be much significant. At 372 GHz, the CSO 10.4-m antenna has a 
half power beamwidth of about 22 . 

Similar setups were used for the para-H30+(l~-2|) line at 307.1924100 GHz (JPL catalogue to 
be found at http://spec.jpl.nasa.gov/), which was observed at CSO in October 2002 and June 2003. 
Tab. lA.l [ present the beam efficiencies that were used for HjO^ data. Pointing was measured every 
two hours and found to be better than 3". 

Observations of the para-D2H+ (li,o-lo,i) line (v = 691.660483(20) GHz; 



Amano & Hirao 



20051) were carried out at CSO between April 2003 and April 2005 under very good weather con- 



ditions (225 GHz zenith opacity less than 0.065). We used the 50 MHz AOS with a spectral resolu- 
tion better than 0.04 km s"' . The observations were performed using the wobbler with a chop throw 
between 150" and 180" according to its stability. The beam efficiency was carefully and regularly 
checked on Mars, Venus, Saturn and Jupiter, and found to be ~ 40%. For more extended sources, 
the beam efficiency has been measured to be ~ 60%. This value has been adopted for 16293E, 
the only source where D2H^ has been detected, assuming that in this case the emission covers an 
area not significantly smaller than the beam. For the other sources, we use rji, = 40% (Tab. |4]i and 
consider a factor of 1 .5 uncertainty in the column density upper Umit value due to the unknown 
source size. Pointing was monitored every 1.5 hours and found to be better than 3 . At 692 GHz, 
the CSO 10.4-m antenna has a half power beamwidth of about 1 1 . 



2.2. Source Selection 



The source list is in Tab. \T\ which reports the coordinates, the Local Standard of Rest velocity 
(Vlsr) at which we centered our spectra, and the distance to the source. 

The selection criteria for starless cores is similar to that described in ICrapsi et al. I (l2005h . 



where sources with bright continuum and N2H^ emission have been selected to include chemi- 
cally evolved cores, where CO is significantly frozen onto dust grains and where HaD^ is thus 
expected to be more abundant. The sample consists of "shallow" cores, with central densities of 
-10^ cm (L 1498, TMC-2, L 1517B, B 68) and more centrally concentrated ones, with central 
densities of -10^ cm"^ (TMC-IC, L 1544, L 183, Oph D, L 429 and L 694-2). 

The star forming regions observed have been selected as being representative of the early phases 
of protostellar evolution, so that any detection of H2D^ will be interesting to compare with star- 
less cores and see if any evolutionary trend can be highlighted. Among the protostellar cores we 
selected: 

(1) NGC 1333 DCO^, clo se to IRAS 4A, in Perseu s, wher e large abundan ces o f deuterated species 



have b een observed (ND3, 



ND2H, 



Roueff" et al. 



2005). 



van der Tak et al. 120021: NH2D, 



Hatchell 



2003 



D2S, 



Vastel et al. 112003 



(2) B 1 is one of the highest colu mn density cores in the Perseus Complex, with active low 



Hirano et al. 



mass star formation going on (e.g. 
as shown by the detection of tripl y deuterated ammo nia (ND3, 
deutera ted hydrogen sulfide (' D2S 



1999b and wit h large deuterium fracti onations. 



Lis et al. 



2002b 



(D2CS. 



Marcelino etal. 112005) 



Vastel et al. 



20061) . doubly 



2003|), and doubly deuterated thioformaldehyde 
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Table 1 : Source sample 



Source 


RA(J2000) 


Dec(J2000) 




Distance 


Name 


[h m s] 


[deg ' "] 


(km s"') 


(pc) 


Starless Cores 


L 1498 


04 10 53.67 


+25 10 18.12 


7.80 


140 


TMC-2 


04 32 44.03 


+24 23 32.56 


6.10 


140 


TMC-IC 


04 41 38.81 


+26 00 21.98 


5.20 


140 


L 1517B 


04 55 18.00 


+30 37 43.84 


5.80 


140 


L 1544 


05 04 17.23 


+25 10 42.70 


7.14 


140 


L 183 


15 54 08.56 


-02 52 48.99 


2.50 


110 


OphD 


16 28 28.56 


-24 19 25.03 


3.40 


165 


B 68 


17 22 38.64 


-23 49 46.03 


3.40 


125 


L429 


18 17 05.53 


-08 13 29.94 


6.82 


200 


L 694-2 


19 41 05.03 


+ 10 57 01.99 


9.60 


250 


Protostellar Cores 


NGC 1333 DCO+ 


03 29 12.10 


+31 13 26.30 


6.00 


350 


B 1 


03 33 20.84 


+31 07 34.27 


6.00 


350 


IRAM 04191 


04 21 56.91 


+ 15 29 46.10 


6.60 


140 


L 1521?" 


04 28 39.80 


+26 51 35.00 


6.48 


140 


OriB9 


05 43 08.17 


-01 15 11.76 


9.10 


450 


NGC 2264G~VLA2 


06 41 11.09 


+09 55 59.01 


8.00 


800 



Embedded low luminosity object recently detected by Spitzer dBourke et al. Il2006h 



(3) IRAM 04191, a very low luminosity Class source in Taurus, driving a powerful outflow, 
but embedded in a dense core which appears to maintain many of the starless core characteristics 
iBelloche & Andre Il2004i . 

(4) L 1521F, initially selected a s a starless core w ith chemical and physical structure similar to 

JCrapsi et al. l2004i . and recently found associated with a L < 0.07 



L 1544 (but different kinematics; 



Lq pr otostellai' object thanks to the high sensitivity of the Spitzer Space Telescope dBourke et al. 



2006h . 



(5) Ori B9, a ma ssive dense core in Orion B, with pecuharly narrow molecular line widths and low 
gas temperature ( Lada et aLlll991 ; Hariu et aLlll993 ; Caselh & Myers |[l995 ), thus an ideal target 
(among massive cores) to detect deute rated species. 



(6) NGC2264-VLA2, studied by, e .g..lWard-Thompson et al 
ing source of the bipolar outflows. 



Girart et al. 



(Il995l) . who found the Class driv- 



( 200C ), who determin ed the gas temperature and 
J2OO2I) . who observed D2CO and 



Loinard et al. 



volume density of the surrounding core, and by 
found an extremely large deuterium fractionation ([D2CO]/[H2CO] = 0.4, equivalent to an enrich- 
ment over the cosmic D/H ratio of more than 9 orders of magnitude). 
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3. Results 



7 



3.1. orfho-H2D+ 



The ortho-H2D^(li,o - li,i) spectra are shown in Fig.[Tl whereas the results of Gaussian fits to 
the hnes are listed in Tab. |2] One striking thing which stands out from the figure and the table 
is the large variation in intensity (factor of 5) and linewidths (factor of 4), and not just between 
starless cores and pr otostellar cores. In column 5 of Tab. |2] we also report the non-thermal line 
width, defined as (see|j 



Myers etal. 111991) 



Avnt = ^A^l^,, - Av2, 



(2) 



where Avobs (= Av in Tab. ^ is the observed line width and Avt is the thermal linewidth of the 
observed molecule, calculated assuming the kinetic temperatures listed in column 3 of Tab. [3] 
It is interesting to note that the ortho-H2D^(li o - li,i) non-th ermal line width is on average 
two times larger than that derived from the N2D^(2-1) data of ICrapsi et al. I (120051) (using the 
same kinetic temperature), the only exception being B 68, where both the N2D^(2-1) and ortho- 
H2D^(li_o - 1 1,0 lines are totally thermally broadened. The larger non-thermal line widths (and the 
moderate optical depths, see Sect.|4|i indicate that the ortho-H2D^ lines are tracing a region within 
the cores with more prominent internal motions than the N2D^(2-1) lin es, even in L 1544, where 
the emission of these two lines has the same morphology and extension ( Vastel et al. 112006a ). 

In the category of "shallow" cores, we have 3 non detections among 4 objects. The only shallow 
core detected in our survey is L 1517B, which is in fact the most compact and centrally concen- 
trated of its class and has one of the narrowest H2D'*' lines observed so far (together with the denser 
core Oph D): 0.4 kms onl y 1 .2 times larger th an the H2D^ thermal linewidth at 10 K, as mea- 
sured with NH3 observations ( Tafalla et al. II2OO4 ). Recent observations carried out with the APEX 
telescope have a ctually revealed a probabl e ortho-H2D^ emission in one of our non-detected shal- 
low cores, B 68 dHogerheiide et al. 1120061) . The detected line is indeed relatively faint (Tmb - 0.2 
K), very narrow (0.3 kms"\ practically thermal), and consistent with our non-detection (Trms - 
0.3 K, see Tab.|2]l. The observations of the other two objects in this group (L 1498, TMC-2) have 
better sensitivities (Tnns ^0.13 K), but they may still hide the H2D^ line if its intensity is similar 
to the one in B 68. 

Among the five most centrally concentrated objects in the starless core sample, four show strong 
(T^^ = 0.7 - 0.9 K) H2D+ emission, whereas TMC-IC has T^b ^ 0.4 K. The line widths span the 
range between 0.4 kms"' (for Oph D) and 0.7 kms"' (for L 429), maybe reflecting contraction 
motions in difl'erent stages of core evolution or large optical depths (although simple analysis seems 
to discard this last hypothesis; see Sect.|4]i. 

Also in the six (young) protostellar cores, the detection rate is quite large, with four H2D^ 
lines detected. L 1521F has a line shape quite similar to that in L 1544 (Av = 0.5 kms"'), but 
the brightness temperature is 1.7 times lower, in agreement with the two times lower deuterium 
fractionation observed (A^(N2D+) /A^(N2 H+) = 0.1 and 0.2, in L 1521F and L 1544, respectively; 



Crapsi et al. 



2004 



Caselli et al. 



2002ah . B 1 and NGC 1333-DCO+ present the largest linewidth 



in the sample, suggesting that the active star formation is probably injecting energy in form of non- 
thermal motions and turbulence. We are confident that the broad line in NGC 1333-DCO^ is not 
a baseline artifact, in particular because both the centroid velocity and the linewidth are coincident 
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Fig. 1: ortho-H2D^(l 10 - li,i) spectra toward the 16 sources of our sample. The units are main 
beam brightness temperature in K (y-axis, assuming a unity filling factor) and velocity in km s"' (x- 
axis). A star in the top right indicates dense cores associated with protostellar objects. The vertical 
dotted line is the Vlsr velocity measured with N2H+(l-0), whereas the vertical dashed Une marks 
Vlsr as measured with the present H2D'*' observations. Within the uncertainties, the two values are 
identical. Note that the strongest emission is present in the densest starless cores (L 1544, L 183, 
L 429, L 694-2) and in the star forming regions L 152 IF, B 1, and NGC 2264G. Note also the large 
variation in linewidth among the various sources. 



(within the errors) with those observ ed in D2S by 
observed by 



Vastel et al. 



van der Tak et al. 



(I2OO3I) (but not with the ND3 line 



2OO2I which remains a puzzle). The H2D+ line in NGC 2264G is 



narrower than in B 1 and NGC 1333-DCO^ and more similar to L 1521F, probably indicating that 
the circumsteUar environment is still quite pristine. 



3.2. para-D2H+ 

The para-D2H^(li o-lo,i) line has been searched for in four sources (B 1, NGC 1333-DCO^, 
NGC 2264G-VLA2, and L 183). Tab. |4] lists the spectral resolution (Avres, column 2), the sys- 
tem temperature (Tsys, column 3) and the integration time (fim, column 4) of the observations. The 
rms noise and the upper limits of the radiation temperature (or brightness temperature, assuming a 
unity filling factor) are in columns 5 and 6, respectively. From these data, we calculated the corre- 
sponding upper limits of the para-D2H^ column density in each source, applying the same method 
as for ortho-H2D^ (see Sect. |4]i and assuming a critical density for the transition of 10'' cm"^, a 
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Table 2: Gaussian fits of the ortho-H2D^(l i o - 1 1,1) lines in the source sample. 



Source 


/ T„,tdV 


Vlsr 


Av 


Avnt" 


Tmb 


Tib 


T 

ims 


Name 


(Kkms"') 


(kms"') 


(kms"') 


(kms"') 


(K) 




(K) 


Starless Cores 


L 1498 












0.6 


0.13 


TMC-2 












0.6 


0.13 


TMC-IC 


0.19±0.02 


5.29+0.03 


0.44+0.06 


0.34+0.09 


0.41 


0.5 


0.12 


L 1517B 


0.14±0.01 


5.85+0.02 


0.40+0.04 


0.29+0.05 


0.33 


0.5 


0.06 


L 1544'' 


0.50±0.03 


7.28+0.02 


0.50+0.03 


0.41+0.03 


0.93 


0.7 


0.11 


L 183 


0.40±0.02 


2.45+0.01 


0.41+0.03 


0.30+0.04 


0.92 


0.7 


0.12 


OphD 


0.28±0.05 


3.40+0.03 


0.39+0.09 


0.27+0.11 


0.67 


0.4 


0.26 


B68 












0.4 


0.32 


B68'^ 


0.078+0.015 


3.36+0.05 


0.33+0.10 


0.0+0.1 


0.22 




0.08 


L429 


0.65±0.05 


6.64+0.03 


0.73+0.06 


0.67+0.07 


0.85 


0.4 


0.21 


L 694-2 


0.50±0.05 


9.56+0.03 


0.53+0.07 


0.45+0.08 


0.88 


0.45 


0.24 


Protostellar Cores 


NGC 1333 DCO+ 


0.25±0.02 


7.65+0.06 


1.49+0.12 


1.44+0.12 


0.16 


0.7 


0.05 


B 1 


0.60±0.03 


6.66+0.03 


1.20+0.08 


1.14+0.08 


0.47 


0.7 


0.09 


IRAM 04191 












0.65 


0.23 


L 1521F 


0.28±0.02 


6.58+0.02 


0.53+0.03 


0.41+0.04 


0.50 


0.6 


0.07 


OriB9 












0.5 


0.25 


Ori B9'' 


0.07±0.03 


9.20+0.06 


0.42+0.15 


0.25+0.25 


0.16 




0.04 


NGC 2264G-VLA2 


0.22+0.02 


4.56+0.03 


0.51+0.04 


0.29+0.07 


0.41 


0.7 


0.11 



Data from 
Data from 
Data from 



Vastel et al. 


1 2006al). 


Hoeerheiide et al. 1 


12006) 


Hariu et al. 


2006 


. 



Une width equal to that measured for t he HaD^ line (Tab. |2l l, except for 16293E (for which we used 
the para-DaH^ line width observed by 



Vastel et al. 



2003b . the kinetic temperature and H2 volume 



densities listed in Tab.[3] and the following parameters for the para-D2H^(l i_o-lo,i) transition: fre- 
quency vo = 691.660483 GHz, Einstein c oefficient for spon t aneou s emission A^i = 4.55x10 s 
rotational constants as given in Tab. 4 of lAmano & Hirao I (l2005h . lower state energy of 50.2 K, 



and degeneracy of the upper and lower levels equal to 9. 

The last column of Tab. |4] shows the ratio between the upper-limit column densities of para- 
DaH^ and the orth o-H2D^ column dens ities listed in Tab. [3] The estimated values are well within 
those calculated bv lFlower et al. I (120041) (see their figure 7) for cloud cores with volume densities 



2x10^ cm ^ and temperature ranges between 10 and 15 K. 
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Table 3: ortho-H2D^column densities. 



Source 


«(H2) 




Ref.'' 


A'(H2)'' 


Ref. 


^ ex 




A'(ortho~H2D+)'^ 


Name 


(cm-3) 


(K) 




(10^2 cm-2) 




(K) 




(10'^ cm-2) 


Starless Cores 


L 1498 


1x10' 


10 


4 


3.2 


7 


7.5,5.5 


<0.3,<0.9 


<0.7,<2 


TMC^2 


3x10' 


10" 


7 


6.0 


7 


8.7,5.9 


<0.2,<0.6 


<0.5,<1.5 


TMC-IC 


1x10* 


7 


8,9 


8.5 


9 


6.8,6.0 


0.4,0.6 


0.9,1.4 


L 1517B 


2x10*^ 


9.5 


4 


3.7 


7 


8.0,5.6 


0.2,0.9 


0.4,1.8 


L 1544 


2x10* 


7 


21 


13 


T 


6.9,6.5 


1.2,1.7 


3.2,4.5 


L 183 


2x10* 


7 


12,22 


14 


12 


6.9,6.5 


1.1,1.6 


2.5,3.5 


OphD 


5x10' 


7 


23 


11 


23 


6.7,5.8 


0.8,2.0 


1.6,4.0 


3 68= 


3x10' 


10 


24,13 


1.4 


7 


8.7,5.9 


0.1,0.3 


0.2,0.5 


L429 


6x10' 


7" 


7 


12 


T 


6.8,6.3 


1.1,1.8 


4.1,6.7 


L 694-2 


9x10' 


7" 


7 


11 


r 


6.9,6.3 


1.1,2.0 


3.2,5.7 


16293Ef 


2x10* 


12 


19 


50 


19 


12,10 


0.5,0.7 


1.2,1.7 


Protostellar Cores 


NGC 1333 DCO+ 


1x10*'' 


13 


1 


31 


14 


12,8.7 


0.03,0.1 


0.35,0.59 


B 1 


3x10* 


12 


2,3 


21 


15 


12,10 


0.1,0.1 


0.9,1.1 


IRAM 04191 


5x10' 


10 


5 


17 


16 


9.2,6.8 


< 0.3, < 0.6 


< 0.8, < 1.5 


L 1521F 


1x10* 


9.3 


6,20 


14 


7 


8.9,7.3 


0.2,0.4 


0.7,1.1 


Ori B9S 


TxlO-* 


10 


10 


4.0 


10 


6.8,4.2 


0.1,0.4 


0.2,0.9 


NGC 2264G-VLA2 


8x10' 


15 


11 


10 


17 


14,9.2 


0.1,0.2 


0.3,0.5 


NGC 1333 IRAS 4A'' 


2x10* 


13 


18 


31 


18 


13,10 


0.03,0.04 


0.2,0.3 



" When more than one reference is listed, the first number refers to the source paper for the volume density 
and the second number refers to the source paper for the kinetic temperature. In the case of L 183, both 
parameters can be found in the two listed papers. 

Values of A'(H2) have been estimated assuming a constant gas temperature of 10 K along the line of sight, 
except for those objects with central temperatures of 7 K (see item "i" below). 

The first value is calculated assuming a critical density «ct=10' cm"^, whereas the second value is for n^r 
= 10* cm ^ 
" Assumed value. 

ortho-H2D^ data from|] 



' ortho-H2D*data from 



Hogerheiide et al. 



Vastel et al. 



1 2006f) . but method described in the text. 



1 2004f) . but method described in the text. 



^ Estimates of Tex, t and A'Cortho - H2D^) based on APEX data from lHariu et al. I ( |2006|) . but using the 

method described in the text. 
^ Data froml! 



Stark et al. 



1999h, but method described in the text. 



' The value of A'(H2) given bv lCrapsi et al. I ( l2005h has been modified to take into account the temperature 
gradient (see Sect. l4.1l for details). 



Refere nces: (1) 
(5) 



Hatchell 



12003); (2) 



Belloche & Andre 



Schnee & Goodman 



Girart et al. 



13" beam; (15) 



Gomez et al. 



i 200C);(12)|Paganietal 



Lis et al. 



1 20051); (9) 



1 2006); (3) 



1 20041) ; ( 6)ICrapsi et al 



Schnee et al. 



Hatchell et al. 



1 19941) . within 20"; (18) 



2004ft ; (13) 



Marcelino et al. 



20041); (7) 



20071); (10) 



12005); (4) 



Crags^^tal 



Hariu et al. 



Hotzel et al. 



1 2005 ), in a 14" beam; (16) 



Tafalla et al. 



1 20041) 



1 20051) , in a 11" beam; (8) 
I993I), in a 40" beam; (11 ) 



l 2002h;(14) 



van der Tak et al. 



Belloche et al. 



(20021), in a 



Stark et al. 



1 2004 ) {N QI2) m easu red in a 13" beam us ing 



1 20021) . in a 1 1" b eam; (17) 



1 199% (mn-,) me asured in a 13" beam) ; (19) 



Lis et al. 



Crapsi et al. 



1 20071) ; (22) 



Pagani et al. 



1 2003); (23) 



Vastel et al. 



2002a data); (20) 



Kirk et al. 



Stamatellos et al. 



1 2003); (24) 



1 20071); (21) 

1 20061) . 



Bergin et al 
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Table 4: p-D2H^ column density upper limits. 



Source 






^int 


T 

-* rms 


Tr 


1h 


f b 
^ ex 


N(p - D2H+) 


picf 


Name 


(km s"') 


(K) 


(min) 


(K) 


(K) 




(K) 


(10" cm-2) 




B 1 


0.083 


3950 


251 


0.068 


<0.20 


04 


11.9 


<0.9 


<1.0 


NGC 1333-DCO+ 


0.083 


3300 


91 


0.068 


<0.20 


04 


12.6 


<0.9 


<2.6 


NGC 2264G-VLA2 


0.021 


1436 


93 


0.045 


<0.14 


04 


144 


<0.2 


<0.7 


L 183 


0.042 


1770 


42 


0.084 


<0.25 


04 


6.9 


<7.9 


<3.2 


16293E'' 


0.10 




103 


0.19 


0.85 


0.60 


10.7 


1.2±0.5 


1.0 


L 1544"= 


0.021 




230 


0.18 


<0.54 


04 


7.6 


<2.5 


<0.8 



^ In radiation temperature units. 

Excitation temperature assuming a critical density of lO' cm"^^ for the observed D2H+ transition. 
" p/o = N(p ara - D^H+yNjort ho - H2D+). 



Data from 
the text. 
" Data from 



Vastel et al. 



1 2004)). The column density has been estimated using the method described in 



Vastel et al. 



2006ah . The column density has been estimated using the method described in 



the text, with the parameters listed in Tab. |3] 



4. Analysis 

4.1. Derivation of the average orfho-HaD^ coiumn densities 

In this section we estimate the average ortho-HiD^ column density in each source, assuming that 
the Une is emitted in homogeneous spheres at the density and temperature respectively quoted in 
the literature and reported in Tab. [3] Values of the H2 column densities are also reported in Tab. [3] 
column 5, and they are used to determine the fractional abundance of ortho- H2D^ (see Sect. 14. 3I) . 
The A^(H2) values for L 1544, L 429 and L 694-2 have been determined bv lcrapsi et al.N2005l) 
assuming constant t emperature. However, a temperature gradienthas been measured toward L 1544 
dCrapsi et al. Il2007h and assumed (because of the sim ilar physical structur e) in L 429 and L 694 



Stamatellos et al. 



2, so t hat A^(H2) needs to be modified. As found by iPagani et al. I (120041) and 
( l2007h . the temperature drop in L 183 and Oph D implies an increase in A^(H2) by a factor of about 
1.4. Given that L 1544, L 429 and L 694-2 have structures similar to L 183 and Oph D, we simply 
multiplied the ICrapsi et al. I (120051) values by the same correction factor (1.4) to account for the 
temperature gradient and this is what is reported in Tab.[3j 

We further assume that the level structure of the ortho-H2D"^ molecule is reduced to a two- 
level system. This is especially justified in starless cores, considering that the first excited state is 
17.4 K above ground and the second one is 110 K. In star forming regions, we assume that the 
observed H2D^ line is arising from gas with characteristics not significantly diflterent from starless 
cores (which is probably true in L 1521F and NGC 2264G-VLA2, see Sect.[3]l. In the case of B 1 
and NGC 1333-DCO^, the broad lines suggest that the embedded young stellar objects have prob- 
ably increased the degree of turbulence in the region and may have locally altered the conditions 
where H2D^ emits. However, we should point out that where the gas temperature increases above 
20 K, and/or where shocks are present, H2D^ should not survive for long, considering that in these 
conditions the backward reaction (1) can quickly proceed and that dust mantles can be either evap- 
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orated or sputtered back into the gas phase, with the consequence of increasing the CO abundance 
and thus the destruction rate of HaD^. 

To estimate the average ortho-H2D"^ column density, we evaluate the excitation temperature Tex 
and the Une optical depth t simultaneously, by solving iteratively the following equations (valid for 

Tbg « E,,/k): 



£ui V «ioi / 



(3) 



1) 



(4) 



-In 



1 



(5) 



Eq. (O is the definition of Tex corrected for the line opacity: in practice, the critical density is 
reduced by the probability that the emitted photon can indeed escape absorption. We ado pted the 
photon escape probability j8 of Eq.|4l valid for an homogeneous sphere ( Osterbrock II 1 989 ). Finally 
the line optical depth t is derived from the observation via equation (|5]l. In the above equations, 
Eai is the energy of the transition {E^xlk^ - 17.4 K, with the Boltzmann constant), «tot is the 
particle volume density (assumed to be 1.2xn(H2), to account for He), Tkin is the gas temperature, 
and JviJesC) and Jy{T\,^) are the equivalent Rayleigh-Jeans excitation and background temperatures. 
The ortho-H2D^ column density is then derived from r: 



Niortho - 



t/ 



Tdv 



(6) 



A key parameter is the critical density of the transition. Recent and unpublished calculations 
by E.Hugo and S.Schlemmer (private communication) suggest a collisional coefficient with ortho 
and para H2 of ~ 10"^ cm^s ' at 10 K, with a very shallow dependence on the temperature. The 
implied critical density is ~ 10^ ^. name ly a factor of 10 lower than assumed in previous work 



Black et al 



199C; 



van der Tak et al. 



2OO5I) . The ortho-H2D^ column densities, calculated both 



with the new and old critical density values, are reported in Tab. [3] The lines are either optically 
thin (e.g. in B 68) or marginally thick (^^1, in L 1544, L 183, L 694-2, L 429). From the table, it is 
clear that a factor of 1 variation in the collisional coefficient implies a change in the column density 
value of factors between 1.4 and 4.5, depending on the volume density and kinetic temperature. In 
the following analysis we use the column densities calculated with the 10^ cm"^ critical density. 

The four objects which show the largest values of Niortho - H2D^) are among the most cen- 
trally concentrated cores in the sample: L 429, L 1544, L 694-2 and L 183. The two other dense 
cores, TMC-IC and Oph D have significantly lower ortho-H2D^ column densities, and this may 
be related to different evolutionary stages. We will further discuss these issues in Sect.|5] 



4.2. Evaluation of the errors 

The estimates of the ortho-H2D^ column densities reported in Tab. |3] suffer from several sources 
of uncertainties. We already mentioned a basic source of uncertainty, that associated with the colli- 
sional coefficient of the transition. In addition to that, the densities and temperatures used to derive 
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the excitation temperatures are also relatively uncertain, not only because of the uncertainty in de- 
riving these values at the centers of the studied sources but also because the H2D^ line emission 
may originate in denser than the quoted average density gas due to the H2D^ abundance distribu- 
tion. In this context, the errors associated with the rms of the observations reported in Tab. |2] are 
certainly the smallest in the error propagation chain. Although it is difficult to exactly quantify the 
error in the determination of the ortho-H2D^ column densities, we estimate here how reasonable 
changes in the gas temperature and density would affect the reported column densities. Increasing 
or decreasing the density by a factor of 2 results in decreasing/increasing the ortho-H2D^ column 
densities by less than 30%. However, a change in the kinetic temperatures of Tab.[3]by 1 K would 
change the ortho-H2D^ column densities by up to a factor of 2 in the coldest objects (because of 
the exponential in the level population equation). 

In summary, considering also the beam efficiency variation between 0.4 and 0.7 (Tab|2]i, the 
ortho-H2D^ column densities reported in Tab.[3]are likely to be uncertain by about a factor of ^2. 

4.3. Correlations 

We have looked for possible correlations between the column density or fractional abundance of 
ortho-H2D^ and physical parameters such as the volume density, the H2 column density, the kinetic 
temperature and the non-thermal line width. No significant correlations have been found, with the 
exception of x(ortho - H2D+) (= N{ortho - H2D+)/A?(H2), with A^(H2) from Tab. O vs. Tkin, for 
which we find (see Fig.|2l): 



Logjc(orf/jo - H2D+) = (-8.6 ±0.3)- (0.16 + 0.03)ridn, 



(7) 



with a correlation coefficient of -0.83. What is causing the ortho-H2D^drop observed between 7 K 
and 15 K? Partially, this may be due to the fact that the warm est sources: (Q may have a smaller 
ortho-H2D^ emitting region than found in L 154 4 (about 60"; IVastel et al. Il2006ah . because of a 



Flower et al. 



20041) and (») are all at distances 



drop of the ortho-H2D^ abundance (see Fig. 6 of |] 
> 300 pc (with the exception of 16293E, which, by the way, follows the trend). Thus, the H2D^ 
emission may be affected by beam dilution (not considered in this study). However, as we will see 
in the next section, variations in the CO depletion factor, linked to the gas density and, in turn, 
to the gas and dust temperature (at least in starless cores, larger temperatures are typically asso- 
ciated with lower gas densities, lower CO depletion factors, and lower deuterium-fractionations), 
can also (at least partially) produce the observed tren d. In any case, a mo re detailed physical and 
chemical structure (such as that recently developed by Aikawa et al.ll2008 ) is needed to investigate 
these points, although the lack of information on the extent and morphology of the H2D^ emission 
prevents us from a more quantitative analysis of the correlations between the gas traced by H2D^ 
and the physical properties of the selected cores. 

In the following we will concentrate on molecular abundances and use both a simple chemical 
model applied to a homogeneous cloud and a slightly more detailed chemical-physical model of 
a centrally concentrated and spherically symmetric cloud to reproduce and interpret the observed 
variations of ortho-H2D^column densities, deuterium fractionation and CO depletion fraction in 
the selected cores. 
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10^ 10® 6 8 10 12 14 16 

n(H2) (cm"') Ti,„, (K) 

Fig. 2: ortho-H2D^ fractional abundances vs. the gas volume density (n(H2), left panel) and gas 
temperature (Tifin, right panel). Empty symbols refer to starless cores, whereas filled symbols refer 
to cores associated with young stellar objects. Note that upside-down triangles are upper limits. 
The line in the right panel is the least square fits to the ji[:(ortho-H2D^) vs. Tkin data, the only 
significant correlation that has been found (see text). 



5. Chemical Discussion 



Roberts et al. 



Aikawa et al. 



(I2OO5I) 



The c hemistry of st a rless c ores and their degree of deuteration has been investigated in de 
tail by 



2003 


), 


Flower et al. 


( 


2004 


), 


Roberts et al. 


(2004 


), 


Walmslev et al. 


Flower et al. 


( 


2005 


), 


Flower et al. 


( 


2006< 


3). 


Flower et al. 


( 


2006b) 



(12004 . 



these models it appears that there are several sources of uncertainty that can profoundly afifect 
the chemistry: (i) surface chemistry (diffusion and reaction rates are still quite uncertain and highly 
dependent upon the poorly known surface of dust grains); (ii) freeze-out and desorption rates 
(binding energies may be changing throughout the core, due to cha nges in grain mantle com - 



Garrod et al. 



2006, 



2007) 



position, and nonthermal desorption processes are poorly known; see 
(iii) the dynamical evolution of dense cores is hard to constrain chemically, and different theo ret 
ical models of star for mation predict significantly different time scales (e.g. 



Shu et al. 



1987 and 



Hartmann et al. 



200 lb ; (iv) the cosmic-ray ionization rate ^ is not well constrained, and our ig- 



norance in the cosmic-ray energy spectrum (especially at low energies) p revents us from rnaking 



Padoan& Scalo 12005 



quantitative estimates on the possible vari ations of within dense cores (see 
for a recent discussion on this point, and iDalgarno II2006I for a more general review); (v) the frac- 
tion of ortho-H2, which affects the deuteration in the gas phase (the back ward reaction (1) pro ceeds 
faster with ort ho-H2, because of its h igher energy compared to para-H2; Gerlich et al. II2OO2 , here- 
after GHR02; I Walmslev et al. II2004I) ; (vi) the dust grain size distribution (if dust grains coagulate 
in the densest regions of st arless cores, the freeze-out rate diminishe s, altering the gas phase chem- 
ical composition; see, e.g. 



Flower et al. 



I2006at 



Vastel et al. 



2006al) ; (vii) the abundance of poly- 



cyclic aromatic hydrocarbons (PAHs), unknown in dark clouds (where observational constraints 



are yet to be found), which may significantly affe ct the electron fraction (ILepp & Dalgarno 



Flower & Pineau des Forets 



2003 



Tielens 



1988 



2005b ; (viii) the fraction of atomic oxygen in the gas 



phase, which is tho ught to be low (mainly to explain the stringent SWAS upper limits on the wa- 
ter abundance; e.g. 



Bergin & Snell 



appear to disprove this dCaux et al 



2002), but observations of dark clouds with the ISO satellite 



1999 



Vastel et al. 



2000; 



Lis et al. 



200 ll) . although the lim- 
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Table 5: The forward rate coefficient is given by a{T lltOfSf. The reverse rate is given by 



Reaction 


Rate 


ff 


P 


7 


Ref. 






■1 1 

cm s 




K 




H+ + HD H2D+ +H2 


k,,k_, 


1.7 10"' 





220 


(1) 


H2D+ + HD ^ D2H+ +H2 


k2, k_2 


8.1 10-'" 





187 


Roberts et al. (2004) 


D2H+ + HD ^ D+ +H2 


kj, k_3 


6.4 10"'" 





234 


Roberts et al. (2004) 


H+ + CO ^ HCO+ + H2 


kco 


1.7 10"' 








RATE06 


H2D+ + CO ^ HCO+ + HD 


kco 


1.7 10"' 








RATE06 


DCO+ + H2 












D2H+ + CO ^ HCO+ + D2 


kco 


1.7 10"' 








RATE06 


DCO+ + HD 












D+ + CO ^ DCO+ + D2 


kco 


1.7 10"' 








RATE06 


H3 + e" — > products 


krffO 


6.7 10"** 


-0.52 





RATE06 


H2D^ + e" — » products 


kj-ff 1 


6.0 10"** 


-0.5 





Sundstrom et al. (1994) 


YioW^ + e" — » products 




6.0 10"*< 


-0.5 





Roberts etal. (2004) 


Dj + e" — > products 


krerf 


2.7 10"** 


-0.5 





Larssonetal. (1997) 



(1) At the low temperatures found in our starless and protostellar cores, this rate corresponds to the Langevin 
limit. 



ited velocity resolution of ISO LWS is a severe limit on these results. In the following section, the 
effects of some of the above parameters on the deuterium fractionation will be presented for the 
simple case of a homogeneous cloud. In Sec. 15.21 a more detailed physical structure and a slightly 
more comprehensive chemical model will be considered to make an attempt on constraining some 
of the unknown parameters. 

5.1. Simple Theory: what affects deuterium fractionation 

Ignoring for the moment the density and temperature structure of molecular cloud cores and any 
gas-dust interaction, which will be considered in Sect. 15.21 we show here simple relations between 



the HaD^/Hj abundance ratio and parameters such as the gas kinetic temperature, the grain size 
distribution, the CO depletion factor and the cosmic-ray ionization rate. We use a simple chemical 
network which includes all the multiply deuterated forms of H3, formed following the reaction 
scheme listed in Tab.|5]and destroyed by CO, electrons (dissociative recombination) and negatively 
charged dust grains (r e comb ination). The adopted rate coefficients are the same as in Tab. 1 of 



Ceccarelli & DominikI (l2005h . with the exception of the reaction of (and deuterated isotopo- 
logues) with CO, and the reaction of H3 and electrons, for which we used the values listed in the 



latest release of the UMIST database (RATE06) available at http://www.udfa.net/ We note that 
the rate coefficient of the H3 -1- CO reaction in the UMIST database is temperature-independent, 
as expected for ion-molecule reactions where the neutral species has a small dipole moment (E. 
Herbst, private communication). 

The steady state equations for this simple system are: 
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a:(H2D+) _ kix(nD)[D2D3 - x(nD)k_2k3] 

x(H+) ~ D1D2D3 - x(HD)(k^3k3Di - k^ikiDj) 

x(D2H+) _ k2x(}iD)D3 

x(H2D+) ~ D2D3 - k^ikixiHD) 

x(D+) _ kixiUD) 

x(D2H^) k-2 + kcoX(CO) + krec3X(e) + kgX(g) ' 

with 



Di - k-i + kcox{CO) + ki-ecix(e) + kgx(g) + k2x(liD) 
D2 — k_2 + kcox(CO) + krec2x(e) + kgx(g) + kjxiHD) 
D3 - A;_3 + ki-ox(CO) + krec3x(e) + kgx(g). 

In the above expressions, ^1, k2, kj, and k^i, k-2, k-3 are the forward and backward rate coefficients 
relative to reactions of , H2D"^, and D2H^, respectively, with HD; x(i) is the fractional abundance 
(w.r.t. H2) o f species a;(HD) = 3 x10'^ is the fractional abundance of HD (assuming that [D]/[H] 
= 1.5x10"^, Oliveira et al. II2OO3 , and that in molecular clouds the deuterium is mainly in the form 
of HD); kco is the rate coefficient of the destruction reaction of all H^ isotopologues with CO; 
;c( CO) = Xcan(CO)/fD, w herc Xcan(CO) - 9.5x10"^ is the canonical abundance of CO as measured 
bv lFrerking et al. I (Il982r) . and fo is the CO depletion factor (1//d is the fraction of CO molecules 
left into the gas phase, see Sect. 15.2.21 ); krec\, kred, krec3 are t he dissociative recombin ation rate 
coefficients for H2 D^, D2H^, and D^, respectively. Following 



Draine & Sutin 



(11987D (see also 



Crapsi et al. II2004I) . the rate coefficient for the recombination onto dust grains has the form: 



-0.5 



^kin 



/t„ = 1.6xl0-^(^^lf^r xfl +3.6x10 
^ \10-**cm/\10K/ \ lOKlO-'^cm 

where flmin (=50 A) is the minimum radius of dust grains in the 



(11) 



Mathis et al. 



(119771) (MRN) size 



distribution (omax = 2.5x 10 ^ cm). Finally, the fractional abundance {x{g)) of dust grains in a MRN 
size distribution is given by: 



jcg = 5.3 X 



ll0-4cm/ \10-**cm/ 



\-2.5 



(12) 



The electron fraction is calculated as in 
reaction scheme of 



Caselli et al. 



Umebavashi & Nakano 



(l2002bl) . using a simplified version of the 



19901) . where we compute a generic abundance of 



molecular ions "mH^" assuming formation due to proton tr ansfer with Y\t and des truction by dis 



Draine & Sutin 



19871) . H+ is formed 



sociative recombination on grain surfaces (using rates from 
as a consequence of cosmic-ray ionization of H2 and destroyed by proton transfer with CO and N2. 
Metals are also taken into account and their fractional abundance has been a ssumed 10"^ (fo llowing 



Lee etal. 111996). 



the initial abundances of "low-metal" models, appropriate for dark clouds; [ 

The deuterium fractionation in species such as HCO^ or N2H^ {Rq = [DCO^]/[HCO^] or 
[N2D^]/[N2H^]) in this chemical scheme is simply given by: 



1 /3x(H2D+) + 2/3x(D2H+) + x(D+) 
x(H+) + 2/3x(H2D+) + 1 /3x(D2H+) ' 



(13) 
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Fig.[3]shows the [H2D+]/[H+], [D2H+]/[H+], [D+]/[H+] andT^o abundance ratios as a function 
of gas temperature, for different values of the depletion factor /d (top panel), the minimum value of 
the dust grain radius in the grain-size distribution, Omin (central panel), and the cosmic-ray ioniza- 
tion rate, ^ (bottom panel). The first thing to note is the sharp drop in the deuterium fractionation 
at temperatures above ^15 K, when reaction ^ and the analogous ones for the formation of D2H^ 
and D3 also start to proceed backward. Between 5 and 10 K, the deuteration ratios increase be- 
cause of the inverse temperature dependence of the H2D^ destruction rate coefficients A;reci, ^rec2, 
kreci and kg (present in the denominator of eqs. [8] 9, and 10). Note the high peaks in the /d - 
50 and 100 [D3]/[H3] curves. This is due to the fact that in regions where most of the neutrals 
are frozen (although here we only talk about CO, we can generalize this statement including in 
the CO-depletion factor all the neutral species which react with the isotopologues), the deu- 
terium fractiona t ion pro ceeds rapidly and H3 is efficiently conve rted into Dt, as already shown by 



Walmslev et al. 



( 120041) in the case of dense cloud cores and by ICeccarelh & DominikI ( l2005h in 
the case of protoplanetary disks. The large abundances of D3 yield large Ru values (> 1). 



Fig.[3]also shows that at temperatures below ^17 K, the deuterium fractionation is very muc h 
dependent upon the CO depletion factor (a well known phenomenon; e.g. lOalgarno & Leppll984t) . 



the gas volume density (see the dotted curves in the top figures), the fractional abundance of dust 
grains and the cosmic-ray ionization rate. In particular, a value of Omin = 5x10"** cm (5 A), coiTe- 
sponds to x(^ ) - 2x1 0~^, which may be regarded as a possible value for the fractional abundance 



of PAHs (e.g. lLepp & Dalgarno 



1988 



Tielens 



20051). Thus, if PAHs are abundant in dense cores 



and they are the main negative charge carriers, the deuterium fractionation is expected to be <0.1, 
because the four H3 isotopologues quickly recombine. The fact that molecules such as N2H^ and 
NH3 show large deuterium fractionations (or large Ro values) in the direction of pre-stellar cores 
and Class sources (Rr, > 0. 1 ; see Fig. |4]l, thus suggests that (negatively charged) PAHs have 
abundances significantly below 10"^. 



The bottom panel of Fig. [3] shows that the larger the cosmic-ray ionization rate the smaller the 
deuteration ratios. This is mainly due to the fact that a larger value of ^ implies a larger electron 
fraction, and a consequently larger dissociative recombination rate (see denominator of eqs. [8] 9, 
and 10). Again, the large deuterium frac tionation observe d in pre-stellar cores and Class objects 
can be used to put upper limits on ^ (see Dalgarno1l2006 ). 



Before proceeding to the next sub-section, we note here that the Rr, values obtained in this 
analysis for typical parameters ((i) fu ^ 10, as typically observed in pre-stellar cores, (ii) Omin - 
50 A, as in the MRN distribution, and (iii) ^ ^ 3x10"'^ s^^) reach ^0.5 for Tki,, < 15 K. This value 
m ay appear too l arge w hen compared to the deuterium-fractionation measured in pre -stellar cores 



bvlCrapsi etal. I ( l2005b . but it is quite close to the [NH2D]/[NH3] ratio found by (l) ICrapsi etal. 



(l2007h in the nuc leus of the pre-stellar core L 1544 using interferometric observations and by (2) 



Pillai et al 



(120071) in Infrared Dark Clouds. However, the results presented in this section apply to 
an ideal homogeneous cloud and are based on "standard" rate coefficients for the proton-deuteron 
exchange reactions H3, H2D^, D2H^-i- HD. In fact, GHR02 have measured slower rates which, if 
adopted, lead to Ru values about a factor of 3 lower compared to those obtained in Fig. [3] This will 
be discussed in the next sub-section. 
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0.5 



0.5 



'r [H2D"]/[H3^]^ 
rf-D= 100 -_ 


J 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ^ 


J 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ^ 


.1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1. 
r Rd ^ 


"1 1 1 1 1 1 1 1 ll 1 1 1 1 1 II 1 ij 

-amin(cm)=5xio-^ - 

5x10^^ 


"1 1 1 1 1 1 1 1 1 1 1 ijiTi 1 


1 1 1 1 1 1 1 1 jTrr^^ 1 1 


"1 1 1 1 1 1 1 1 1 1 1 1 MlT^ 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

(s"^) = 3xl0-'^ - 
3xl0"'^\ '. 




TTT 1 1 1 1 1 1 1 1 1 [Trri 


TTTT|Ti III iTTfiX 1 

^-u-JIl:^t^, , 1 


1 1 1 1 1 1 1 1 1 1 1 1 imTT 
" iTT^ 



5 10 15 20 25 10 15 20 25 10 15 20 25 10 15 20 25 
Tkin (K) Tkin (K) T,i„ (K) Tkin (K) 

Fig. 3: [H2D+]/[H+], [D2H+]/[H+], [D+]/[H+] and Rd abundance ratios as a function of the gas 
temperature (Tkin) in a dense cloud with uniform volume density n(H2) = 10^ cm"^. (Top row) 
The abundance ratios are plotted against Tkin for different values of the depletion factor /d (=1, 
5, 10, 50 and 100), with fixed values of Omin = 50 A and ( - 3x10 s"'. The dotted curves 
are for /d = 10 clouds with n(H2) = 1x10^ cm"-' (top dotted curve) and n(H2) - 1x10'* cm"^ 
(bottom dotted curve). Note the large increase of the [DJ ]/[H3 ] ratio for /d = 50 and 100. (Central 
panel) Abundance ratios vs. r^n for different values of a^i„ and values of /d and ( fixed at 10 
and 3x10"'^ s"', respectively. (Bottom panel) Abundance ratios vs. Tkin for different values of the 
cosmic-ray ionization rate, ^. Here, /d = 10 and a^in = 50 A. 



5.2. Simple Chemical-Physical Model 



In this section, our estimates of ortho-H2D^ column densities are correlated with the deuterium 
fractionation and the depletion factor, previously measured in the same objects, and simple chem- 
ical models are used to investigat e the observed variati ons among the vario us sources. The mode l 

(l2006ah . and first applied by|( 



d2002bh 



used is similar to that described inlVastel et al. 



Caselli et al. 



in the case of L 1544, but with the deuterium fractionation chemistry and rate coefficients as de- 
scribed in the previous sub-section. We consider a spherical cloud, with a given density and tem- 
perature profile, where dust and gas are pr esent. Initially (besid es H2), CO and N2 are present in 
the gas phase with abundances of 9.5x10"^ jFrerking et al. 1 19821) and 3.75x10"^, respectively. The 
abundance of N2 assumes that 50% of the nitrogen is in atomic form (but no Nitrogen chemistry 
is considered, except for the N2 adsorption/desorption onto/from dust grains and the formation 
and destruction of N2H^ and N2D^). This is a totally arbitrary choice, but the N2 abund ance is 



extremely uncertain (see e.g. IStepnik et al. 



2003 



Flower et al. 



2006b: 



Maret et al. 



2006) and in 
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any case, its variation in the gas phase only affects, in our simple model, the absolute abundance 
of N2H^, without significantly affecting the N2D"^/N2H"^ column density ratio. Atomic oxygen has 
not been included in the chemistry, in order to avoid one extra (uncertain) parameter in the model. 
It is worth to point out here that adding atomic oxygen to the chemical network lowers the deu- 
terium fractionation , if its binding energy (also not well constrained) is sufficiently low (see e.g. 



Caselli et al. Il2002bl and their discussion in Sect. 3.2). 

The dust grain distribution follows MRN and we assume a gas-to-dust mass ratio of 100. 
However, the size of the minimum dust radius, a^i„, has been increased by an order of magni- 



-""■-'5 ^^^^^ ^. 

Flower et al. 


2005 


Bergin et al. II2OO6 


Flower et al. 


2006b 


Vastel et al. 


2006a) 



adopted here (5x10 ^ cm) lowers the freeze-out rate by a factor of 5 compared to the MRN value 
(by changing th e surface area of dust grains), and it is th e "best-fit" value for the L 1544 chem- 



ical model (see 
(^fi-eeze(O) is givcn by: 

^freeze (0 — 



Caselh et al. 



2002b 



Vastel et al. 



2006al) . The freeze-out time scale of species ; 



5 I < vth(0 > nil 
2 X lO^r 



(^min \ / ^/jxiax \ 

10-5cm/ "llO-5cm/ 



-0.5 



A(0' 



o.5/10K\° 



lO^cm- 



(14) 



where S is the sticking coefficient (^ 1, as recently found by iBisschop et al. Il2006r) . < Vth(0 > is 
the average thermal velocity of species /, «h is the total number density of H nuclei («h = n(H) + 
2n(H2)), A{i) is the atomic mass number of species /, and 



(4.88 X 10-^^ + 4.66 x 10"^^) x (a^'^^ 



"max /' 



(15) 



is the grain surface area per H nucleon in the MRN distribution (see also IWeingartner & Draine 
200 lb . The electron fraction is calculated as in Sect. 15. II 



The binding energies for CO and N2 have been taken from lOberg et al. I (120051) . assuming that 
the mantle composition is a mixture of CO and H2O ice (£'d(CO)/^b = 1 100 K and £'d(N2) = 0.9 
£'d(CO)). The cosmic-ray ionization rate used here is f = 1.3x10 s but we have also changed 
it to explore the effects on the chemistry (see next subsections). Different density structures have 
be en considered (see b elow) and the (gas - dusQ ) temperature profile is similar to the one found 
by I Young et al. 1120041 and parametrized so that: 



Tdust ~ rga.s - 3 X [8.7 - log(n(H2)] K. 



(16) 



The minimum (maximum) allowed temperature is 4 K (14 K). We also consider models with tem- 
perature profiles increasing inwards, to simulate the heating of the embedded young stellar object. 
In this case, we assume a central temperature of 50 K at a distance of 80 AU, and a radial depen- 
dence r for r > 80 AU. If the temperature drops below the one described in eq.(fT6b, the latter 
value is used. 



This a ssumption of simi lar dust and gas temperatures is only valid if the densities are larger than ~10^^ 



Goldsmith 



20011) . However, in the range of temperatures typical of low-mass dense cores, this 



cm ' (e.g. 

approximation does not significantly change the results of our model. 
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CO and N2 can freeze-out (with rates given by l/ffreeze, see eq. [T4b and return to the gas 



Hasegawa & Herbst 



pha se via thermal desorption o r cosmic-ray impulsive heating (foUovying iHasegawa et al. 
and 



1992 



19931). The abundance of the molecular ions (N2H^, HCO^, Hi and their 



deuterated isotopologues) are calculated in terms of the instantaneous abundances of neutral species 
(as sumption based on th e short time scale of ion-chemistry, compared to the depletion time scale; 
see Caselli et al. IboOZb for details). 



5.2.1 . Niortho - H2D+) vs. the observed 



In Fig. m the column density of ortho-H2D^ is plotted as a function of the observed deuterium 
fractionation ratio {Rd). Rd is equivalent to Af(N2D^)/A^(N2H ^), in the case of star less cores (plus 
L 1521F), and this value has been taken from the survey of ICrapsi et al. I ( l2005h . In star form- 
ing regions, the N2D"^/N2H'^ column density ratio is not available, and other colum n density ra- 



tios have been used: (i) A?(NH2D)/A?(NH3) for NGC 1333-DCO+ dHatchell 



2003) and for B 1 



(IRoueffetal. KOOSh : (ii) ^/N(D^CO)/N(ihCO) for NGC 2264G VLA2 dLoinard et al. hooA No 
deuterium fractionation estimates are avai lable for IRAM 0419 1 . Given that NH^ an d N2H"^ appear 
to trace similar zones of dense cores (e.g. 



Benson et al. 



1998 



Caselli et al. 



2002c ), and that both 



derive from the same parent species (N2), one expects that the D-fractionation observed in the two 
species is also similar (and linked to the theoretical Rd in eq.[T3]). However, it is not obvious that 
formaldehyde is actually tracing the same region (indeed H2CO is centrally depleted in the two 
starless cores studied by Tafalla et al. 2006, unlike N2H^ and N2D^). In fact, Fig.|4]shows that the 
deuterium-fractionation in NGC 2264G is the largest one in the whole sample, probably suggest- 
ing that different deuteration mechanisms (beside the H3 fractionation) may be at work for H2CO. 
One possibility is that surface chemistry is needed to ex plain the observed am o unt of deuterated 
forma ldehyde and methanol, as originally discu ssed by 



(119981) . and more recently by 



Parise et al. 



Charnlev et al 



(Il997h . 



Ceccarelli et al. 



(12006D. 

In the left panel, each data point is labelled with the corresponding name, whereas the right 
panel shows the same data points with model curves superposed (see below). The first thing to note 
is that, on average, dense protostellar cores (filled symbols) have lower N{ortho - H2D^) values 
than starless cores, but on average they show quite large deuterium fractionations (especially in 
the case of NGC 2264G VLA2, where R^ is coming from measurements of doubly deuterated 
formaldehyde, as already mentioned). Another thing to note is that there is not any clear correlation 
between N{ortho - H2D^) and the observed /?d- To investigate this unexpected result, we used the 
m odel described in Se c. 15721 and simulate an evolutionary sequence, similar to what has been done 
in lCrapsi et al. I (l2005h in their Fig. 5. 



We consider Bonnor-Ebert (BE) spheres with density structures analogous to the radial (cylin- 
dr ical) density profile of t he contracting disk-like cloud at different stages of evolution in the model 
of ICiolek & Basu I (I2OOOI) . namely those at times t - ti (=2.27 Myr and central densities nci(H2) 
= 4x10"* cmr\ ti (=2.60 Myr, and nc2(H2) = 4x10^ cm^^), fs (=2.66 Myr, and naini) = 4x10*' 
cm-\ U (=2.68 Myr, and ticAmi) = 4x10^ cm"^), and (=2.684 Myr, andnrf(H2) = 4x10^ cm"^) 



The B E density profile is reasonably well reproduced by the parametric formula dTafalla et al 

«.(H2) 



2002) 



n(r) = 



(1 + (r/ro)«) 



(17) 
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Observed Rr, Observed R, 



Fig. 4: ortho-H2D^ column density as a function of the observed deuterium-fractionation Rd (see 
text for details). Filled circles are protostellar cores, whereas empty circles are starless cores. 
Downward arrows denote upper limits in our estimate of the column density. The names of the 
cores associated with each mark in the figure are shown in the left panel. The right panel shows 
the same plot, where theoretical predictions from chemical models (see text) are superposed. Solid 
curves are for models which use the "standard" rate coefficients for the proton-deuteron exchange 
reactions, whereas dashed curves are for models adopting the smaller GHR02 rate coefficients. 
Each curve has six points (filled (empty) diamonds for the models using the "standard" (GHR02) 
rate coefficients), which correspond to model cores with different central density (from 4x10^* cm 
to 4x10** cm"^, from left to right). Different curves are for different ortho/para HiD^ ratios, from 
0.03 to 2.0. The cosmic-ray ionization rate has been fixed to 1.3x10"'^ s"'. 



where ro = 13,000, 3,000, 800, 300, and 80 AU for fi, fi, fs, f4, and fs, respectively, and a = 2. 
The five chemical models (i.e. the CO and N2 depletion chemistry in the five model clouds with 
different physical structure) have been run for a time interval given by (f, - fi), for / = 2, 3, 4, 
and 5, whereas the fi model has been run for 2.27 Myr. We also consider a model cloud with the 
same density structure as the fj model, but with a temperature profile resembling that of a centrally 
heated protostellar core (model t^a), with a central temperature of 50 K and a temperature gradient 
proportional to r (as mentioned above). The abundance profiles of HiD^, N2H^ and N2D^ 
calculated by the models have been convolved with 22", 26", and 17" FWHM antenna beams, 
respectively, to simulate the observations and calculate the column densities toward the core center 
(as well as off peak). 

The results of these models are the small diamonds in each of the curves in the right panel of 
Fig. m with f 1 lying on the left-end and fsa on the right-end of the curve. Solid curves represent 
models with standard rate coefficients for the proton-deuteron exchange reactions, whereas dashed 
curves models use the about 3 times smaller GHR02 rates (see previous sub-section). The different 
curves correspond to models with different values of the o/p ratio ofHjD^ (o/p-H2D"^), from 0.03 
(bottom curves) to 2.0 (top curves). As discussed by 



Flower et al. 



( I2OO4I) . the o/p ratio is a sensitive 
function of the o/p H2 ratio and, ultimately, of the gas temperature (see their Fig. 6) and at Tg^s < 
~15 K, it changes fr om ^0.03 to values probably larger than one (this last statement is valid if 
the curve in Fig. 6 of iFlower et al. I (120041) is simply extrapolated at temperatures lower than 9 K, 
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the minimum value in the figurj^). In all curves, the models show a slight decrease in both the 
H2D'*" column density and in the deuterium fractionation. In fact, at such high central densities: (i) 
H2D^ is efficiently converted into D2H^ and D:^ , thus lowering the total H2D^ column along the 
line of sight; (ii) N2 significantly freezes onto dust grains, so that the N2D^/N2H^column density 
ratio - our measure of the observed - traces regions away from the center, where the density is 
lower and the temperature is higher. If a protostar is present, the H2D^ column density increases 
again because of the less efficient transformation of into D^, while Ru decreases (see also 
Fig. [3]) because of the less abundant D3. From the comparison between the data and the model, 
one is tempted to conclude that indeed variations in the o/p-H2D^ ratio (and ultimately in the gas 
temperature of the central few thousand AU of dense cores) can explain the observed scatter among 
the cores. It is interesting to note that the two pre-stellar cores with high values of Ru and relatively 
low ortho-H2D^ column densities are both embedded in the Ophiuchus Molecular Cloud Complex: 
Oph D and 16293E. Chemical abundances observed in these two cores are consistent with a lower 
H2D^ o/p ratio, suggesting possible (slightly) larger kinetic temperatures. 

Fig.|5]shows two other attempts to interpret the data. The left panel considers exactly the same 
mo dels as in Fig. [4] but now the o/p-H2D^ratio is fixed at 0.5 (consistent with dense gas at 9 K; 



Flower et al. 



20041) . whereas the free parameter is the cosmic-ray ionization rate which is 



varied from 6x10 s ' (bottom dashed and solid curves) and 3x10 " s ' (top curves). We note 



that v ariations in the cosmic-ray ionization rate are known to exist in the Galaxy dvan der Tak et al 



20061). The four pre-stellar cores with the largest N{ortho - H2D+) values (L 492, L 1544, L 694-2, 



and L 183) can all be reproduced by f2-f4 (t\-t2) models with the GHR02 (standard) rates and f - 
1x10"''' s"' (with L 429 (L 183) being the most (least) dynamically evolved). Significantly lower 
values of ^ (< 6x 10"* s ') appear to be required in the protostellar and Ophiuchus cores. However, 
an alternative way to reproduce these data, without requiring extremely low ^ values, is to lower the 
o/p-H2D^ ratio, as found before (the thick dashed curves represent models with o/p-H2D^ - 0.3 
and ^ = 6x10"* s '). In B 68, only the low o/p-H2D^ model curve (at early evolutionary times) 
can reproduce the low ortho - H2D^ column densities and deuterium fractionations, the lowest in 
the sample. 

In the right panel of Fig. |5] we consider five clouds with structure as in the f, {i- 2, 5) 
models (see eqn. [TtI i and for each of them we follow the chemical evolution, checking the results 
after 10^, lO"', 10^, and 10^ yr. This allows us to explore how different chemical ages can change 
the gas composition and avoid the problem of fixing the (unknown) chemical evolution time as in 
the right panel of Fig. |4] As for the left panel, the o/p-H2D^ ratio has been fixed at 0.5, except 
for the thick curves, representing the fja protostellar cloud models, where o/p = 0.1, assuming that 
the gas has been (slightly) heated compared to the pre-stellar cores. The cosmic-ray ionization 
rate is fixed at 1.3x10 s The four ortho-H2D^- rich pre-stellar cores can be reproduced by f3 
models, with (chemical) ages between lO'* and 10* yr, when GHR02 rate coefficients are used. On 
the other hand, the protostellar cores and the two Ophiuchus cores are better matched by the more 
dynamically evolved (centrally concentrated) fga and fg models, respectively, after about 10^-10^ 
years of (chemical) evolution. 



^ At a density of 2x10* cm ^, temperature Tgas = 8 K, grain size 0.1 yum, the 0/P-H2D+ is about 1.2 (Pineau 
des Forets, priv. comm.) 
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Fig. 5: (Left panel) Same models as in the right panel of Fig. |4]but with a fixed value o/p-H2D^ 
(=0.5), and various cosmic-ray ionization rates (from 6 to 30 xlO s '). The thick dashed curve 
connects model results using GHR02 rates, ^ - 6x10 s ' and o/p-H2D^= 0.3. This is an attempt 
to reproduce the low ortho-H2D^ column density values observed toward the /?D-rich protostellar 
cores and the two Ophiuchus pre-stellar cores Oph D an d 16293E. (Righ t pane l) Plot showing 



Ciolek & Basu 



(I2OOOI) cloud at times 



predictions of models with density structures equal to the 
ti (i - 2, 3, 4 and 5) and coiTesponding central densities nd = n, (see text) at different stages 
of chemical evolution (at t = 10^, 10"^, 10^, and 10* yr, from bottom left to top right). The o/p- 
HaD^ ratio has been fixed at 0.5. The thick curves refer to the model of a protostellar envelope at 
difi^erent evolutionary times, where o/p-H2D^= 0.1, to simulate a possible (slight) increase of the 
gas temperature due to the central heating source. 



5.2.2. N(ortho - H2D+) vs. the observed /d 



From the models described in the previous subsection, one can directly derive the abundance of 
CO within each cloud model as a function of cloud radius and, as before, N(CO) is obtained, after 
integrating the CO number density along the line of sight and smoothing the results with a 22" 
beam width (to simulate observations carried out at the IRAM 30m antenna, where most of the 
cores have been observed). From the model column density, the CO depletion factor, /d(CO) is 
easily calculated using the expression: 

:can(CO) 



/d(CO) = 



(18) 



N(CO)/N(U2y 

where .Xcan(CO) is t he "canonical" or "u ndepleted" abundance of CO (assumed here equal to 
9.5 xl0~^, following iFrerking et al. Ill982l but known to be uncertain by a factor of about 2, see 



e-g- 



Lacy et al. II 1994 . 



To compare our model predictions with the data, we collect from the literature the values of ob- 
served /d(C O) and plotted them versus N(ortho - H2D^) in Fig. |6| The majori ty of the /d(CO) dat a 



comes from 



Crapsi et aL 



(l2005l). except fo r: (i) TMC-IC (/d(CO) = 6.9, from 



Schnee et al. 12007) 



(ii)L 183 (/d = 5, from lPagani et al. l2005l): Cm) NGC 13 33-DCO+ (/d = 12, from 



2OO2I): (\v^ B 1 rfn(CO) = 3.2, from 



Belloche & Andre 



Caselli & Myers 



Lis et al. 



J0rgensen et al. 



2004]); (vi) OriB 9 (/d(CO) = 3.6, from iHariu et al. 



2002bl) : (v) IRAM 04191 (fn( CO) = 3.5 from 

for A^(H2), and 



2006 



1995L for N(CO)). Also, for L 1544, L 429, L 694-2 and Oph D, the new val- 



ues of A^(H2), adopted to take into account the temperature structure (see Tab. [3), imply different 
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Fig. 6: N(ortho - HaD^) vs the observed CO depletion factor, /d(CO). The same models f, (i=l, 
5) described for the previous figures are used here, varying the cosmic-ray ionization rate ( (left 
panel) and the o/p-HaD^ ratio {right panel) within the range of Fig.|4]and|5] The thick dashed curve 
refer to (GHR02) models with 0/P-H2D+ ratio of 0.1 and ^ = 6x10 s As found before, varia- 
tions in the o/p-H2D^ ratio, plus differences in the physical structure, provide a way to reproduce 
the whole observed spread in the data. The labels in the left panel correspond to the source names: 
A=L 1498, B=TMC-2, C=TMC-1C, D=L 1517B, E=L 1544, F=L 183, G=Oph D, H=B 68, 
I=L 429, J=L 694-2, KsNGC 1333-DCO+, L=B 1, MsIRAM 04191, N=L 1521F 0=Ori B9. 



values of /d (larger by a factor of about 1 .4, the ratio between th e new and old A^(H2) values, as 
explained in Sect. 14. It from what reported bv lCrapsi et al. I (120051) . The data and model results are 



shown in Fig. |6] In general, the presence of embedded young stellar objects appears to lower the 
H2D'^ column density, without affecting much the amount of CO freeze-out, which is probably still 
large in the high density and cold protostellar envelopes. The possible (small) temperature increase 
caused by the central heating, is thus not sufficient to significantly release CO back in the gas 
phase, but it can affect the o/p-H2D^ ratio (and, consequently, the ortho-H2D^ column density), as 
discussed in the previous sub-section. 

In the left panel of Fig.|6] the data are compared with the same models described in Fig.|5](left 
panel; ti, i=l, 2, ...5), where the cosmic-ray ionization rate is varied from 6x10 s ' to 3x10 
s~'. The o/p-H2D^ ratio has been fixed at 0.5, except for the thick dashed curve, where o/p = 0.1, 
value probably more appropriate for protostellar cores (see previous sections). Now, f2-f3 models 
(depending on the rate coefficient values adopted for the proton-deuteron exchange reactions), with 
( ^ 10 s can explain the observations toward the ortho-H2D^- rich pre-stellar cores, with the 
exception of L 183, where the low fu value suggests a younger dynamical phase, consistent with 
what found in the previous sub-section (we also note that L 183 and TMC-IC appear to have similar 
ages, but with different cosmic ray ionization rate and/or o/p-H2D^ ratio, see right panel of Fig.|6|. 
Lower o/p-H2D^ ratios are needed to reproduce the protostellar and Ophiuchus cores, as found for 
the N(ortho - H2D+) vs. Ro relation. 

The right panel of Fig. |6] shows the same set of data and models, but now f has been fixed at 
1.3x10 s \ whereas the o/p-H2D^ ratio has been changed as in Fig. ID As already noted, L 1544, 
L 429, L 694-2, and L 183 data are best reproduced by large values of the o/p-H2D^ ratio, and the 
appropriate physical structure is that of ti-t^ models , similarly to the left panel (with L 183 being 
the least evolved). This implies cores slightly more evolved than what found in Fig.|4](right panel) 
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and Fig. |3] (left panel), where (standard rate) models between fi and f2 were preferred. The small 
discrepancy can be understood if the predicted CO depletion factor is too large compared with 
observations. Reasons for this could be: (i) we are missing an important desorption mechanism 
(besides the cosmic-ray impulsive heating and the thermal desorption, the latter not being efficient 
at the temperatures of these objects); (ii) unlike our spherical and isolated model cores, real cloud 
cores are embedded in molecular clouds where CO is quite abundant. Thus, the observed "extra" 
gaseous CO may be part o f the undepleted ma terial accreting onto the core from the surrounding 



molecular cloud (see also 
TMC-IC). 



Schnee et al. 



2007 



for a similar conclusion in the particular case of 



6. Conclusions 

Low-mass dense cloud cores have been observed with the CSO antenna at the frequency of the 
ortho-H2D^(l 1 0-1 1,1) line. The main conclusions of this work are: 

1) In starless cores, the line has been detected in 7 (out of 10) objects. The brightest lines (Tmb - 
0.7-0.9 K) are observed toward the densest and more dynamically evolved starless cores (L 1544, 
L 183, Oph D, L 429, and L 694-2), where the deuterium fractionation and the CO depletion factor 
are largest. Non-detections are found in cores less centrally concentrated and dense than the rest of 
the sample (L 1498, TMC-2, and B 68). In B 68, recent APEX observations have detected a faint 
ortho-HiD^ line with intensities consistent with our upper Umit. 

2) Significant differences in line widths are observed also among starless cores, with the naiTowest 
lines (Av ^ 0.4 km s ') found in TMC-IC, L 1517B, Oph D, and L 183, whereas L 1544, L 429, 
and L 694-2 show Av = 0.5-0.7 kms"'. 

3) The ortho-H2D^(l 1 - 1 1,1) line has been detected in 4 out of 6 protostellar cores, where we find 

< 0.5 K (the largest value observed is toward L 152 IF, which hosts a very low luminosity object 
recently detected by Spitzer). The broadest lines are observed toward the two cores in Perseus 
NGC 1333-DCO+ and B 1, where Av ^1 kms The ortho-H2D+(li,o - li,i) Hnes have broader 
non-thermal widths than N2D^(2-1) lines, even in L 1544, where the two transitions have similar 
extension and morphology. 

4) The ortho-H2D^ column densities range between 2 and 40x10'^ cm"^ in starless cores and 
between 2 and 9 xlO'^ cm"^ in protostellar cores. Thus, protostars in the earliest stages of their 
evolution appear to have already changed the chemical structure of their envelopes, by lowering the 
ortho-H2D"^ fractional abundance but not their deuterium fractionation. This is probably due to (a 
few degree) variation of the gas temperature, which strongly affects the o/p-H2D^ ratio. A similar 
effect is also found in the two Ophiuchus cores, suggesting a (slightly) larger gas temperature than 
in the other (mainly Taurus and Perseus) cores. 

5) Simple models suggest that variations in the gas kinetic temperature, CO depletion factor, grain 
size distribution, cosmic-ray ionization rate and volume densities can largely affect the H2D^ abun- 
dance relative to H3 . In particular, gas temperatures above 15 K, low CO depletion factors and large 
abundance of negatively charged small dust grains or PAHs drastically reduce the deuterium frac- 
tionation to values inconsistent with those observed toward pre-stellar and protostellar cores. 

6) Plots of the ortho-H2D^ column density vs. (i) the deuterium fractionation observed in species 
such as N2H^ and NH3, and (ii) the observed CO depletion factor, show a large scatter. The data 



26 P. Caselli et al.: Survey of ortho-H2D*(l i,o-l 1,1 ) in dense cloud cores 

can be reproduced by chemical models of centrally concentrated cores assuming variations in the 
o/p ratio of H2D'*" (ultimately linked to kinetic temperature variations). Changes in the cosmic- 
ray ionization rate between 6 and 30x10 '**s ' can also explain part of the scatter, but not those 
objects with large deuterium fractionations and low ortho-H2D^ column densities (such as the 
Ophiuchus pre-stellar cores and the protostellar ones), for which a decrease of the o/p-H2D^ ratio 
is required. The most deuterated and H2D^-rich objects are better reproduced by chemical models 
of centrally concentrated cores with ric = a few X 10^ cm"^ and chemical ages between lO'* and 
10^ yr. Protostellar cores, plus the two Ophiuchus cores, are better matched by lower o/p-H2D^ 
ratios, more dynamically evolved (central densities >10^ cm"^) models and chemical ages of ^10^- 
10"* years. To better constrain dynamical and chemical ages, the rate coefficient for the proton- 
deuteron exchange reaction needs to be well defined. 

7) The para-H30^(lj^-2|) upper limits are consistent with radiative transfer calculations if the 
fractional abundance of H^O^ is $10"^. The para-D2H^(l i,o-lo,i) upper limits provide an upper 
limit of the para-D2H^ to ortho-H2D^ column density ratio (= p/o). We find p/o<l in all the 
sources (except for L 183 and NGC 1333-DCO^, where p/o < 3), consistent with ch emical model 



predi ctions of high density (2x10* cm ^) and low temperature (r^n < 10 K) clouds dFlower et al. 



2004) 



More accurate determinations of temperature and density profiles, as well as observations of 
the para-H2D^(lo,i-lo,o) line at 1.37 THz, are sorely needed to put more stringent constraints on 
gas-grain chemical processes in dense cloud cores. 

Appendix A: The search of para-HsO^ 

The H^O^ line was observed to investigate the chemistry of oxygen in dense cores and, with the 
help of chemical models, to put some constraints on the oxygen abundance, which, analogously to 
CO, significantly affects the deuterium fractionation. 

We searched for para-H30^(l^-2p in seven dense cores but only upper limits were measured. 
Tab. lA.ll lists the results of this search, including the rms noise and the corresponding upper limits 
of the column density, which have been calculated in two different ways: (i) using the volume 
density and kinetic temperature values listed in Tab.[3](A^i, column 6), and (ii) assuming a volume 
densi ty of 10^ cm~^ and Tj^m = 10 K (A^2, column 8). In both cases, the RADE^il code has been 
used dvan der Tak et al. The observed upper limits (see Tab. lA.ll ) are compatible with para- 



HjO^ fractional abundance upper limit s of 10 ^ according to calculation made using the Ratran 



code dHogerheiide & van der Takll2000h 
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